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The preparation of monolayers of graphite (graphene) 
has led to exciting discoveries associated with the unique 
electronic structure of this two-dimensional system. fP, ^ 
Transport experiments are usually performed in a field- 
effect device geometry with a graphene fiake separated 
from a gate electrode by a dielectric spacer. [3, 4 Doped 
Si is commonly used as the gate material and Si02 as the 
dielectric. Substrate inhomogeneities [SJ |6] and trapped 
charges [71 E] lead to local variations of the electrostatic 
potential at the Si02 surface which results in uncon- 
trolled local electron- and hole-doping of the graphene 
sheet, so-called electron-hole puddles. ^|8^ 

With a layered graphene-like honeycomb structure and 
a lattice constant only 1.7% larger, insulating [10] h- 
BN would be the ideal substrate for graphene based 
devices; [ll] the interaction between /i-BN and graphene 
sheets is so weak that even if their lattices are forced to 
be commensurate, the characteristic electronic structure 
of graphene is barely altered. For an incommensurate 
stacking of graphene on /i-BN, the perturbation is even 
smaller, leading to an extremely well-ordered graphene 
layer with a very high carrier mobility. [12] (TS] [14] 

Like graphene, /i-BN layers can be prepared by 
mechanical exfoliation. [15 Cleaved layers can also be 
thinned to a single layer with a high-energy electron 
beam. [16] Alternatively, /i-BN layers can be grown by 
chemical vapor deposition (CVD) on transition metals 
such as Cu or Ni, using precursors such as borazine 
(B3N3H6) or ammonia borane (NH3-BH3).[17l [18] With 
a proper choice of growth conditions, homogeneous ul- 
trathin /i-BN only 1-5 atomic layers thick can be grown. 
Moreover, graphene can be grown by CVD on top of a 
/i-BN layer adsorbed on a metal substrate, [19 which is 
ideal for field-effect devices. 

A field effect is created by applying a voltage differ- 
ence between graphene and the metal substrate, result- 
ing in charge accumulation or depletion in the graphene 
layer. We study this so-called electrostatic doping as 
a function of the thickness of the /i-BN layer and the 
applied voltage using first-principles (DFT) calculations 
for a Cu(lll)|/i-BN|graphene structure. In the absence 
of an applied voltage, a spontaneous charge transfer 
across the /i-BN layer occurs between the metal sub- 
strate and the graphene, leading to an intrinsic doping 
of graphene. This transfer, which is driven by the differ- 
ence in graphene and metal work functions, is strongly 



modified by the charge displacements resulting from the 
weak chemical interactions at the metal |/i-BN and at the 
/i-BN I graphene interfaces. 

By varying the applied voltage, the doping level can 
be controlled. Experimentally the position of the Fermi 
level in graphene has a square-root like dependence on 
the voltage. [20, 8,|2ll|ll[l4] We wih demonstrate that 
this behavior is reproduced by our DFT calculations on 
Cu(lll)|/i-BN|graphene and develop an analytical model 
that describes how the Fermi level depends on the ap- 
plied voltage and the /i-BN layer thickness. 

Computational Details. We use DFT at the level 
of the local density approximation (LDA), within the 
framework of the plane-wave PAW pseudopotential 
method, [22 as implemented in VASP. [27, "2T, "2F In our 
previous work we found that the LDA gives a reasonable 
description of the geometries resulting from the weak 
interactions between /i-BN layers and between a /i-BN 
and a graphene layer. [JLl The change in the electronic 
structure of such systems in an external field is also de- 
scribed wen.[2ll[27l[28] The LDA also describes the inter- 
action between graphene and metal (111) surfaces very 
reasonably. [29l [30l [3TJ [32l [33] We expect it to provide a 
description of the interaction between /i-BN and metal 
(111) surfaces of similar quality. In contrast, we found 
that PW91 or PBE GGA functionals (incorrectly) pre- 
dict essentially no interlayer binding between /i-BN or 
graphene layers, or between metal (111) surfaces and h- 
BN or graphene sheets. 

The Cu(lll)|/i-BN|graphene structures are modeled in 
a supercell periodic in the z direction with six Cu atomic 
layers in the (111) surface orientation, a slab 1-6 layers 
thick of /i-BN, a graphene monolayer, and a vacuum re- 
gion of ~ 15 A. A dipole correction is applied to avoid 
spurious interactions between periodic images of the slab. 
We choose the lattice constant of graphene equal to its 
optimized LDA value a = 2.445 A, and scale the in- 
plane lattice constants of Cu(lll) and /i-BN so that the 
structure can be represented in a 1 x 1 graphene surface 
unit cell. [11] [29l [30l |31^ The positions of the atoms in 
graphene and /i-BN, and of the top two atomic layers of 
the metal surface are allowed to relax during geometry 
optimization. The effect of the strain in the Cu and h- 
BN structures on the electronic structure is very small. 
For instance, the work function of Cu(lll) is increased 
by a mere 0.04 eV and the density of states at the Fermi 
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Figure 1: The Fermi level shift A£^f in graphene versus 
the thickness D of the /i-BN film. The symbols represent 
DFT results for Cu(lll)|/i-BN|graphene structures with 
2-6 layers of /i-BN. The dashed line represents the model 
of eq. [3) 

level is unaltered. An electric field applied across the 
slab is modeled with a sawtooth potential. [34 

In the most stable configuration of a /i-BN monolayer 
on Cu(lll), the nitrogen atoms are adsorbed on top 
of Cu(lll) surface atoms and the boron atoms occupy 
hollow sites. The calculated equilibrium separation be- 
tween /i-BN and the Cu surface is 3.11 A. The graphene 
and /z-BN layers are stacked as in Ref. 11 with a calcu- 
lated graphene- /i-BN equilibrium separation of 3.21 A, 
and a separation between /i-BN layers of 3.24 A. The 
equilibrium structures are changed negligibly by the ap- 
plied electric field, at least for field strengths up to 0.5 
V/A. The interaction of graphene with /i-BN is so weak 
that the characteristic linear band structure of graphene 
about the conical points is essentially preserved; a small 
band gap of ~ 50 meV is induced in graphene if the /i-BN 
and the graphene lattices are commensurate. If they are 
incommensurate, the gap seems to disappear. p!3 | fT2 | IT4] 
In the following this point will not be important. 

Intrinsic Doping. We monitor the doping of graphene 
in the Cu(lll)|/i-BN|graphene structure by calculating 
the change of the Fermi level AE-p with respect to neutral 
graphene. A computationally convenient way [30l |3T] of 
characterizing AE-p is in terms of the difference between 
the work functions W of the metal | dielectric | graphene 
stack and Wq for free-standing graphene. 



AEp = W-Wg, 



(1) 



Negative (positive) values of A£^f then correspond to re- 
type (p-type) doping. Figure [l] shows AEp as a function 
of the number of /i-BN layers. These numbers are calcu- 
lated without an external field which means that in the 
Cu(lll)|/i-BN|graphene structure the graphene sheet is 
intrinsically doped with electrons. AEp decreases with 
the number of /i-BN layers. We will model this thickness 
dependence below. 



Figure 2: Plane-integrated electron density difference 
An{z) for a Cu(lll)|/i-BN| graphene structure with five 
/i-BN layers. The positions of the top layer of Cu atoms, 
the bottom and top /i-BN layers, and of the graphene 
layer are indicated by vertical dashed lines. Each of the 
three An{z) curves represents a different value of the 
external electric field £^ext- 



The intrinsic doping originates from electrons that are 
transferred from the Cu electrode to the graphene sheet 
across the /i-BN layer. However, the calculated work 
functions of Cu(lll) and graphene are I^m = 5.25 eV 
and Wq = 4.48 eV, respectively. If charge transfer is 
driven by the difference between the Cu substrate and 
graphene work functions alone, then establishing a com- 
mon Fermi level would require transferring electrons from 
graphene to Cu and the result would be p-type doping, 
i.e. a positive value of AE-p^ at variance with the results 
shown in figure [l] The /i-BN layer must therefore play a 
non- trivial role. 

The interaction between two materials at an inter- 
face generally results in the formation of an interface 
dipole. The latter can be visualized in terms of the 
electron rearrangement at the interface characterized by 
the electron density of the entire system minus the elec- 
tron densities of the two separate materials (using iden- 
tical, frozen atomic structures). As only the depen- 
dence perpendicular to the interface is relevant, it is 
convenient to work with plane-integrated electron den- 
sities n{z) = ff n{x^y, z) dxdy where the integration is 
over the surface unit cell. The electron displacement 
in a Cu(lll)|/i-BN|graphene stack is then described by 
An{z) = nM|BN|Gr(^) - ^m(^) - riBNiz) - nGr(^). The 
result for a Cu(lll)|/i-BN| graphene structure with five 
layers of /i-BN is shown in figure [2] 

Dipoles are clearly visible at the interfaces between 
the different materials. The largest interface dipole is 
between Cu and /i-BN. Electrons are piled up on the Cu 
surface implying depletion close to the /i-BN surface. A 
similar effect is observed in the physisorption of organic 
molecules on metal surfaces, [35l |36] or even in the ad- 
sorption of noble gas atoms on metal surfaces. [37J There 
it is attributed to Pauli exchange repulsion between the 
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adsorbate and the substrate, which results in a "push- 
back" of electrons towards the "softer" material, in this 
case the metal substrate. [37] This mechanism for inter- 
face dipole formation is quite general, and we speculate 
that the dipole at the Cu — /i-BN interface has a simi- 
lar origin. Between /i-BN and graphene, one observes a 
similar but significantly smaller interface dipole. These 
dipoles are truly localized at the interfaces; their sizes 
do not depend on the number of /i-BN layers as long as 
there is more than one. Both dipoles can also be ob- 
tained in separate calculations for the two interfaces, i.e. 
one for /i-BN adsorbed on Cu(lll), and one for graphene 
adsorbed on /i-BN. 

An interface dipole layer results in a discontinuity in 
the potential energy perpendicular to the interface. The 
potential energy step at an A|B interface can be obtained 
from an A|B slab calculation as the difference between 
the work functions on the A and the B sides of the slab. 
In practice, a single /i-BN layer on top of the Cu(lll) 
surface is sufficient to calculate the potential energy step 
Acu|BN at the Cu(lll)|/i-BN interface. Similarly, the 
potential energy step Abn|c ^it the /i-BN| graphene in- 
terface can be determined with a system comprising two 
/i-BN layers and the graphene sheet. 

We find Acu|bn = 1.12 eV, and Abn|c = 0.14 
eV. Without additional charge transfer the difference 
between the Fermi levels in Cu and graphene in the 
Cu(lll)|/i-BN|graphene structure would be 

Vq = Wm-Wq- Acu|bn - Abn|c (2) 

which we calculate to be Vq = —0.49 eV. To achieve equi- 
librium requires transferring electrons between Cu and 
graphene to set up an electrostatic potential that com- 
pensates for Vq. The sign of Vb indicates that electrons 
are transferred from Cu to graphene, which results in re- 
type doping, i.e. a negative value for AE-p^ in agreement 
with figure [T] 

External Eield. The charge transfer leads to an intrin- 
sic electric field across the /i-BN slab, which polarizes the 
/i-BN layers. This polarization is clearly identifiable in 
figure [2] as small oscillations of An{z) in the /i-BN region 
in the absence of an external electric field (red line). It 
can be eliminated by applying an external electric field 
that opposes the intrinsic field. With £^ext = +0-1 V/A 
the total internal electric field is zero, and the charge dis- 
tribution in the /i-BN slab becomes identical to that of a 
free-standing /i-BN slab (black line). The blue line in fig- 
ure |2] shows An{z) resulting from reversing the external 
field ^ext = -0.1 V/A, that increases the polarization of 
the /i-BN slab. 

An external field can also be used to control the 
position of the Fermi level, i.e. the concentration of 
charge carriers, in graphene. [3l |4] Calculated Fermi 
level shifts for a range of /i-BN slab thicknesses and ex- 
ternal electric field values are plotted in figure |3] The 
curves for different /i-BN slab thicknesses have a similar, 
highly nonlinear, shape. Similar shapes of the position 
of the Fermi level as a function of a gate voltage have 
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Figure 3: Fermi level shift as a function of an external 
electric field for Cu(lll)|/i-BN| graphene structures with 
2-5 layers of /i-BN. The symbols indicate calculated DFT 
values; the lines represent the model of eq. [3j 




Figure 4: Schematic drawing of the metal|/i- 
BN I graphene structure and the electrostatic potential 
(energy) across the structure. £^int represents the to- 
tal electric field across the /i-BN slab of thickness d; the 
other symbols are explained in the text. 

been observed in scanning tunneling spectroscopy (STS) 
experiments, [ 20\ [H [131 [H] as well as in work function 
measurements. [21] 

The points on the curves where AEp = correspond 
to the charge neutrality level of graphene, i.e. to un- 
doped graphene. At these points the external field ^ext 
is equivalent to a potential difference Vg across the /i-BN 
slab that compensates for Vq. We expect |£^ext| \Vg\/d^ 
where d ^ D the thickness of the /i-BN layer. The exter- 
nal field strength corresponding to the charge neutrality 
level should then decrease monotonically with increasing 
slab thickness, as is indeed observed in figure [3| 

Model. To understand the intrinsic doping (fig. [T]), 
as well as the external field effect (fig. [3| quantita- 
tively, we develop the following analytical model whose 
parameters are shown in figure [4| The work func- 
tion of the metal I dielectric I graphene stack is given by 
W = Wm — AV^ with AV the total potential difference 
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across the stack. We model this potential difference as 
AV{d) = -eEi^td + Acu|bn + Abn|c, where £;int is the 
total electric field inside the dielectric, and e > is the 
elementary charge. This field can be related to an exter- 
nally applied electric field by eo(£^ext — ^^int) = with 

the dielectric constant of the dielectric layer and a the 
surface charge density of the graphene sheet. These re- 
lations can be used in eq. [l] to derive a first expression 
for A£^F in terms of a. 

A second relation between AE-p and a is obtained 
by noting that charge in graphene is introduced by 
(de)populating states away from the charge neutrality 
point, a = e J^^^ D{E)dE^ and that the density of 
states near the conical points is well described by a lin- 
ear function, D{E) = Do\E\/A, with Dq = 0.09/(eV2 
unit cell) [30 and A = 5.18 the area of a graphene 
unit ceh. This then gives a = sign{AEp)AE^ eDo/{2A). 
Combining the two relations between A£^f and a gives 



A^F = ± 



y/l^2aDod\V^-Vo\/i^-l 
aD^d/n 



(3) 



where Vg = — e£^ext<^/^7 Vb is given by eq. [2j and 
a = e'^/eoA = 34.93 eV/A. The sign of AE^ is deter- 
mined by the sign ofVo — Vg. Equation [3] describes how 
the Fermi level in graphene depends on the gate voltage 
Vg (or the external electric field ^ext) and the thickness 
d of the dielectric layer. We define d as the separation D 
between the top Cu layer and the graphene sheet minus 
a correction do, because most of the displaced charge is 
localized between Cu and /i-BN and between /i-BN and 
graphene, see figure |2j We have used the value do = 2.4 
A. [31 In principle the dielectric constant of the /i-BN 
layer is a (weak) function of its thickness. We have 
used the constant value = 2.72, which we calculate for 
a two-layer /i-BN slab by putting the isolated slab in an 
external electric field [34 . We then determine the inter- 
nal field by calculating the macroscopic average V(z)[3S^ 
of the electrostatic potential and differentiating it. The 
dielectric constant is then defined as the ratio between 
the internal and external field. 

The results of this model are given in figures [l] and [3j 
The agreement with the results from the DFT calcula- 
tions for the Cu( 111) I /i-BN I graphene structures is very 
good. Note that the model has no adjustable parame- 
ters, as quantities such as Vq (eq. |2| are obtained from 
separate calculations. The model correctly describes the 
intrinsic doping and its dependence on the /i-BN layer 
thickness, fig. [l] as well as the dependence of the doping 
on the external field, fig. [3] 

Summary. The doping of graphene in Cu(lll)|/i- 
BN I graphene structures is studied by monitoring the 
Fermi level shift by means of first-principles DFT calcu- 
lations. We predict that graphene is intrinsically n-type 
doped for sufficiently thin /i-BN layers. This is due to 
the potential difference between Cu and graphene aris- 
ing from substantial dipole layers at the Cu|BN and the 
BN I graphene interfaces, as well as from the work func- 
tion difference between Cu and graphene, eq. [2] The 



doping level decreases with increasing /i-BN layer thick- 
ness, and approaches zero for thick layers. It can be 
varied by applying an external electric field and the re- 
sulting shift of the Fermi level has a modified square-root 
like dependence on the field. For thick dielectric spac- 
ers, and in the absence of work- function-difference and 
interface-dipole terms, this is similar to what is observed 
in experiments. [201 El [H [H [21] For very thin dielec- 
tric layers, these interface terms are predicted to play 
an important role and both the /i-BN layer thickness de- 
pendence as well as the field dependence of the doping 
can be described quantitatively by an analytical model, 
eq. [3] The parameters of this model can be determined 
experimentally or obtained from DFT calculations on in- 
dividual surfaces (the work functions of the metal and 
of graphene), on interfaces (the interface dipoles formed 
at the Cu(lll)|/i-BN and the /i-BN| graphene interfaces), 
or on slabs (the dielectric constant of a /i-BN layer). 
Graphene field-effect devices using a thin layer /i-BN gate 
dielectric should exhibit an intrinsic doping, and a depen- 
dence of the Fermi level on the applied gate voltage that 
is described by eq. [3] 
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